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Abstract. The 500 central pc of the Galaxy (hereafter GC) exhibit a widespread gas component with a kinetic temperature 
of 100-200 K. The bulk of this gas is not associated to the well-known thermal radio continuum or far infrared sources like 
Sgr A or Sgr B. How this gas is heated has been a longstanding problem. With the aim of studying the thermal balance of 
the neutral gas and dust in the GC, we have observed 18 molecular clouds located at projected distances far from thermal 
continuum sources with the Infrared Space Observatory (ISO). In this paper we present observations of several fine structure 
lines ([Oi] 63 and 146 yum, [Cn] 158 yum, [Si n] 35 yum, [S i] 25 yum and [Fe n] 26 yum), which are the main coolants of the gas 
with kinetic temperatures of several hundred K. We also present the full continuum spectra of the dust between 40 and 190 yum. 
All the clouds exhibit a cold dust component with a temperature of ~ 15 K. A warmer dust component is also required to fit the 
spectra. The temperature of this dust component changes between 27 and 42 K from source to source. We have compared the 
gas and the dust emission with the predictions from J-type and C-type shocks and photodissociation region (PDRs) models. We 
conclude that the dust and the fine structure lines observations are best explained by a PDR with a density of 10 3 cirT 3 and an 
incident far-ultraviolet field 10 3 times higher than the local interstellar radiation field. The fine structure line emission arises in 
PDRs in the interface between a diffuse ionized gas component and the dense molecular clouds. The [Cn] 158 yum and [Si n] 
35 yum lines also have an important contribution from the ionized gas component. PDRs can naturally explain the discrepancy 
between the gas and the dust temperatures. However, these PDRs can only account for 10-30% of the total H 2 column density 
with a temperature of ~ 150 K. We discuss other possible heating mechanisms for the rest the warm molecular gas, such as 
non-stationary PDRs, X-ray Dominated Regions (XDRs) or the dissipation of supersonic turbulence. 
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1. Introduction 

The central 600 pc of the Galaxy (hereafter the Galactic center, 
GC) contain up to 10% of the neutral gas of the Milky Way. In 
spite of the high gas surface density in the GC (several orders 
of magnitude higher than in the disk of the Galaxy) the star 
formation rate is 10 times lower than that in the disk and the 
star formation efficiency is only similar to that in the disk. This 
is probably due to the particular physical conditions of the GC 
molecular clouds, which are dense (~ 10 4 crrT 3 ) but have tur- 
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bulent linewidths of ~ 30 km s and are relatively warm. As 
shown by the NH3 observations of Hiittemeister et al. (1993), 
the GC clouds present a wide range of temperatures from -15 
K to ~ 150 K. They also showed that the warm temperatures are 
found all over the GC, not only in the vicinity of Sgr A and Sgr 
B. Regarding the dust, the large scale emission has been stud- 
ied in the photometric far-infrared (between 40 and 250 /jm) 
survey by Odenwald and Fazio (1984). Their map is relatively 
similar to the radio continuum maps, showing some discrete 
sources over an extended emission. They derived a dust tem- 
perature of < 30 K. Recently, Pierce-Price et al. (2000) have 
derived an homogeneous dust temperature of 20 K in the GC 
from 450 and 850 /im data. Clearly, the dust temperature in the 
GC molecular clouds is considerably lower than the gas tem- 
perature. The thermal structure of the GC clouds is difficult to 
explain but it is required a mechanism that heats selectively the 
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gas and acts over large regions. Radiative heating mechanisms 
are usually ruled out with the argument that most of the gas 
does not seem to be associated with thermal continuum sources 
and that the dust temperature is much lower than that of the gas. 
Several alternative heating mechanisms have been discussed in 
the literature as cosmic-rays or magnetic heating (Morris et al. 
1983, Gusten et al. 1981, Huttemeister et al. 1993). However, 
the most promising heating mechanism is shocks or the dissi- 
pation of supersonic turbulence (Wilson et al. 1982, Gusten et 
al. 1985). 

We have selected a sample of 18 sources located in various 
parts of the GC with the aim of having further insight in the 
physical properties of the GC interstellar medium. The sources 
were selected as molecular peaks in CS and SiO large scale 
maps (Bally et al. 1987, Martin-Pintado et al. 1997, 2000). 
They are located far from thermal radio continuum or far- 
infrared sources in order to study the thermal balance of the 
interstellar gas in the sources where the origin of the high tem- 
peratures is more difficult to explain. The typical sizes of the 
clouds are 5-10 pc. Rodriguez-Fernandez et al. (2001a) have 
measured H2 densities of (5-20) 10 3 cm 3 and column den- 
sities of ~ (1 — 6) 10 22 cirT 2 in these sources. The mass of 
the clouds ranges between 10 4 and 10 5 M Q . Therefore, they 
are representative of the GC molecular clouds (see Miyazaki 
& Tsuboi 2000 for a statistical study of the GC clouds popula- 
tion). 

We have observed these clouds with the LWS and SWS in- 
struments onboard the Infrared Space Observatory (ISO) since 
they are perfectly suited to study the thermal balance of the 
GC clouds. For instance, H2 pure-rotational lines observations 
have allowed to measure the total amount of warm gas in the 
GC clouds (Rodriguez-Fernandez et al. 2000, 2001a). On aver- 
age, the column density of warm H2 is ~ 30% of the H2 col- 
umn density traced by CO. However, in some sources it can 
be much higher (up to 100 %). The warm H2 is mainly gas 
with a temperature of 150 K. On the other hand, fine structure 
lines of different ions allow us to probe the low density ion- 
ized gas and to estimate the properties of the ionizing radiation 
(Rodriguez-Fernandez and Martin-Pintado 2004). The whole 
GC region is permeated by radiation with an effective tempera- 
ture of ~ 35000 K. ISO also gives us the possibility of studying 
the dust continuum spectrum around the emission peak (be- 
tween 40 and 190 /mi) and fine structure lines as [Oi] 63 /mi, 
[Cn] 158 /mi or [Si 11] 34 /mi, which are the main coolants of 
the neutral gas with kinetic temperatures of a few hundred K. 
This paper is organized as follows: Section 2 presents the ob- 
servations and explains the data reduction procedures. In Sects. 
3 and 4 we report on the results regarding the dust continuum 
and the fine structure lines. Section 5 revisits the H2 emission 
and Sect. 6 is devoted to a discussion of the implications of 
those results. We summarize our conclusions in Sect. 7. 

2. Observations and data reduction 

We have taken LWS continuous spectra (43-197 //m) in grat- 
ing mode (LWS01). The spectral resolution of this mode is 
0.29 /mi for the first five detectors (43-93 /mi) and 0.6 /mi 
for the rest (84-197 /an). The aperture of the LWS is ~ 80" x 



[Fe ll]25 prn [Si ll]34 /zm [Si l]25 /tm 



M + 0.94-0.36 




-1000 1000 -1000 1000 -1000 1000 



M+0.83-0.10 




-1000 1000 -1000 1000 -1000 1000 



M + 0.58-0.13 




-1000 1000 -1000 1000 -1000 1000 




-1000 1000 -1000 1000 -1000 1000 

E 
u 




-1000 1000 -1000 1000 -1000 1000 



M + 0.21-0.12 




-1000 1000 -1000 1000 -1000 1000 



M+0.1 6-0.10 




-1000 1000 -1000 1000 -1000 1000 



M-0.15-0.07 




-1000 1000 -1000 1000 -1000 1000 



M-0.32-0.19 




-1000 1000 -1000 1000 -1000 1000 



M-0.55-0.05 




-1000 1000 -1000 1000 



velocity (km/s) 
Fig. 3. Sample of the lines observed in the SWS02 mode 

80 . At the resolution of the LWS01 mode a number of fine 
structure lines are detected. Those arising exclusively from 
Hn regions are presented in Rodriguez-Fernandez and Martin- 
Pintado (2004). Here we present the 63 and 145 /mi lines of 
Oi and the 158 //m line of Cn. We have also taken [Si 1] 25.2 
/mi, [Fe 11] 25.9 /mi and [Si n] 34.8 /mi spectra with the Short 
Wavelength Spectrometer in grating mode (SWS02 mode). The 
line wavelengths as well as the telescope aperture and the ve- 
locity resolution at those wavelengths are shown in Tabled In 
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Fig. 1. Position of the observed sources (black diamonds) on the 20 fim image by the MSX satellite (grey levels) and the C O 
map by Dahmen et al. (1997, contours). M+3.06+0.34 and M+2.99-0.06 (located in the Clump 2 complex) are not shown in this 
figure 

Table 1. Integrated flux densities of the lines observed in the SWS02 and LWS01 modes. Numbers in parentheses are the errors 
of the last significant digit 



Line 


SI 


Fell 


Sill 


Oi 


Oi 


Cii 


AQim) 


25.24 


25.98 


34.81 


63 yum 


145 yum 


158 /im 


Aper. ( x ") 


14x27 


14x27 


20x33 


86 x 86 


70 x 70 


68 x 68 


Res. ( km s -1 ) 


230 


230 


207 


1400 


1240 


1140 


Units W cirr 2 


i n-20 

10 


1 n-20 
10 


i n-20 

10 


1 n-19 
10 


i n-19 

10 


i n-19 
10 


M-0.96+0.13 


<0.79 


<0.55 


38(2) 


24.0(11) 


<2.4 


65(3) 


M-0.55-0.05 


<0.90 


<1.6 


158(5) 


77(4) 


<14.0 


242 (7) 


M-0.50-0.03 


<0.59 


1.8(4) 


116(5) 


60(2) 


<8.0 


164(7) 


M-0.42+0.01 


<0.56 


1.8(2) 


148(4) 


45(3) 


<7.2 


145(6) 


M-0.32-0.19 


<0.80 


1.0(4) 


106(4) 


59.9(9) 


<7.2 


142(6) 


M-0. 15-0.07 


<0.64 


3.9(3) 


196(5) 








M+0.16-0.10 


<0.65 


<1.1 


121(3) 


39(3) 


<11.0 


114(8) 


M+0.21-0.12 


<0.68 


13.1(6) 


580(20) 


83(6) 


<10.0 


183(6) 


M+0.24+0.02 


<0.62 


1.3(3) 


146(4) 


53(2) 


<18.0 


140(7) 


M+0.35-0.06 


<1.0 


1.7(4) 


154(3) 


67.1(9) 


<8.6 


168(8) 


M+0.48+0.03 


<0.89 


1.0(2) 


33.6(10) 


31.6(9) 


<9.9 


107(7) 


M+0.58-0.13 


<0.77 


<0.91 


42(2) 


27.8(11) 


<9.4 


100(5) 


M+0.76-0.05 


<0.58 


<0.51 


23.9(10) 


24.9(13) 


<8.5 


51(7) 


M+O.83-0.10 


<1.2 


<0.66 


20(2) 


27.4(13) 


<4.7 


61(4) 


M+0.94-0.36 


<0.6 


<0.67 


5.5(13) 








M+ 1.56-0.30 


<0.71 


<0.49 


5(2) 


10.1(6) 


<2.0 


25.1(6) 


M+2.99-0.06 


<0.89 


<0.51 


17.0(8) 


8.9(6) 


<1.7 


36(2) 


M+3.06+0.34 


<0.70 


<0.49 


2.9(6) 


9.6(8) 


<1.5 


18.8(13) 



addition, the [Oi] 63 and 145 fim and the [Cn] 158 fim lines 
have also been observed with the LWS Fabry-Perot (LWS04 
mode). The velocity resolution in this mode is ~ 30 km s _I and 
the wavelength calibration uncertainties expressed in units of 
velocity are 5 km s _1 for the [Oi] 63 and [Cn] 158 fim and 15 
km s _1 for the [Oi] 145 fim. 

All of the observations were processed using the off-line 
processing (OLP) software version 10. Further reduction was 
done with ISAP. The data reduction consists basically in drop- 
ping bad data points, shifting the different scans taken with 
each detector to a common level and averaging all the data for 
each detector. In addition, with ISAP we have defringed the 
LWS01 spectra to get rid of the well-known interference pat- 
tern. In the overlapping regions of the LWS detectors we have 
determined that the discrepancies in the fluxes measured by 
the different detectors are lower than 15%. We have not tried 
to shift the different detectors to obtain a smooth continuum 
spectrum since there is no reliable way of doing this. The abso- 
lute flux calibration uncertainties of the LWS01 mode are lower 
than ~ 20%. The dust continuum spectra are shown in Fig|2] 

In Figs. and @] we show a sample of the lines observed 
in the SWS02 and LWS01 modes. Table □ contains the line 
fluxes as derived by fitting Gaussian curves to the lines. The 
absolute flux calibration uncertainties of the SWS02 observing 
mode are less than 30 %. Figure[5]shows a sample of the spectra 
obtained in the LWS04 mode. The results of fitting Gaussian 



curves to the lines (fluxes, linewidths and the velocities of the 
center of the lines) are listed in Table[2] The line flux calibration 
uncertainties of this mode are ~ 20%. At the resolution of the 
LWS04 mode all the lines have been resolved. Most of these 
have very broad lines, with linewidths up to ~ 150 km s . The 
central velocities and the profiles of the different lines detected 
in a given source are in agreement with each other. The LWS04 
and LWS01 lines fluxes are in agreement within the calibration 
uncertainties. 

3. Dust 

We have fitted the dust continuum spectra using the Levenberg- 
Marquardt method. It is not possible to obtain a good fit to the 
spectra with just one gray body. Thus, we have used a model 
with two gray bodies in which, 1^, the total intensity radiated at 
wavelength A is defined as: 

h = (1 - e- T " ) B A {T V ) Q t , + (1 - e-^O B A (T C ) fi c (1) 

where subindex c and w refer to a cold and a warm compo- 
nent, respectively. B,\(T) is the Planck function at temperature 
T, Q. is the solid angle of the emission, and the dust opacity 
at wavelength A. Since the dust emission in the GC is extended 
we have used Q w = Q c = Q.lws (Qlws represents the LWS 
beam). To obtain the dust opacity at a given A, we have extrap- 
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Table 2. Integrated flux densities (in units of 10 18 Wcm 2 ), widths and peak velocities (in km s ) of the lines observed in the 
LWS04 mode. Numbers in parentheses are the errors of the last significant digit 



Source 


Cn 157.74 yum 




Oi 63.18 fim 




Oi 145.52 /im 






Flux 


Av 


Vhel 


Flux 


Av 


Vhel 


Flux 


Av 


Vhel 


M-0.96+0.13 


1.9(2) 
2.6 


61 
150 


o 

-100 


<1.4 


(100) 




<0.3 


(100) 




M-0 55-0 05 


5.0(7) 


83 


-100 


2.4(6) 


100 


-110 


<0.5 


(200) 






8.2(6) 


120 


20 


1.1(6) 


71 


26 








M-0 50+0 03 


12.5(1000) 


110 


-140 


4.1(150) 


110 


-160 


<2 


(200) 






7.2(1000) 


85 


-28 


0.8(140) 


40 


-48 








M-0.42-0.01 


5.1 
11.2 


87 
135 


-146 
-30 


3.4(7) 


150 


-130 


<0.4 


(200) 




M-0 32-0 1 9 


5.7(6) 


60 


o 


1.0 


60 


20 


<1.3 


(110) 






10 


200 


-65 


2.8 


170 


-65 








M-0. 15-0.07 




















M+0.16-0.10 


10.2(8) 


122 


20 


3.2(4) 


130 


39 


<0.6 


(150) 




M+0.21-0.12 


20.8(9) 


117 


16.9 


8.0(10) 


93 


36 


0.65(14) 


42 


20 


M+0. 24+0.02 


18.3(13) 


106 


14 


5.1(6) 


108 


33 


2.1(4) 


104 


20 


M+0 35-0 05 


19.9(14) 


87 


27 


7.1(9) 


75 


36 


1.1(2) 


73 


36 


M+0.48+0.03 


11.2(5) 


88 


30 


2.0(4) 


50 


38 


0.8(2) 


95 


43 


M+0. 58-0. 13 


12.8(7) 


120 


39 


4.9(5) 


109 


51 


0.55(7) 


61 


67 


M+0.76-0.05 


2.5(200) 
3.6(200) 


59 
95 


-14 
68 


<1.8 


(140) 




<0.5 


(140) 




M+O.83-0.10 


4(4) 
3.9(4) 


78 
72 


-8 
90 


<2 


(140) 




<0.6 


(140) 




M+0.94-0.36 




















M+ 1.56-0.30 


<1.2 


(150) 




<0.8 


(150) 




<0.3 


(150) 




M+2.99-0.06 


2.4(2) 


74 


-14 


<0.8 


(100) 




<0.3 


(100) 




M+3.06+0.34 


<0.8 


(100) 




<0.8 


(100) 




<0.6 


(100) 





olated the 30 /mi opacity (T30) using a power law with exponent 

t a = r 30 (2) 

In principle ft varies from 1 for amorphous to 2 for crys- 
talline silicates. Using sub-millimeter data, Pierce-Price et al. 
(2000) and Lis et al. (2001) have found higher /3 values (2.5- 
3.4) in some GC sources. These authors interpret their finding 
as a consequence of the presence of icy mantles in the dust 
grains. We have fitted the dust temperatures (T w and T c ) and 
the opacities (t3o„, and T3 Uc ) of the two dust components for 
different values of ft in the range of 1-3. Tableland Fig. [6] 
show the results of fitting the spectrum towards M+0.21-0.12 
(source with the hottest dust) and M+0.76-0.05 (source with 
the coldest dust). Fits of similar quality (as determined from 
X 1 ) can be obtained for ft values in the range of 1-3. The best 
fits have been obtained with /3 = 1 for both the cold and the 
warm components. Fits with higher /3's present x 1 values just 
-15% higher. 

The fit results for th )k and t^) C are strongly dependent on 
f3 c and j3 w . However, the total opacity is dominated by the cold 
dust component (t w /t c ~ 1/50; Table QJ- On the other hand, 
one should note that fitting the data with high /3s requires very 
high dust opacities (T30 ~ 100) that imply a visual extinc- 
tion, A v , of ~ 1000 (r 30 = 0.0UA v ; Draine 1989). These 
huge extinctions are orders of magnitude higher than those esti- 
mated from any other tracer (Rodriguez-Fernandez and Martin- 



Pintado 2004 have found values of Ay=30-60). Therefore, in 
the observed wavelength range one can rule out /3 values equal 
or higher than 2. This result is in agreement with that obtained 
by Goicoechea et al. (2004) for the Sgr B2 envelope. The high 
/? values found with sub-millimeter observations may be due to 
an intrinsic dependence of f3 with the wavelength (Bernard et 
al. 1999). 

Table |5] also shows that in contrast to T30, T w and T c are 
almost independent of f} c and /3 W and can be determined within 
2 or 3 K. Therefore, Table|4]only contains the dust temperatures 
of the cold and warm components in all the sources for /fe in 
the range of 1-3. The temperature of the cold dust component is 
14-18 K for all the sources. In contrast, the temperature of the 
warm component varies from 23-25 K for M+0.76-0.05 to 36- 
42 K for M+0.21-0. 12. We conclude that the GC clouds exhibit 
two dust components, one with a temperature of 14-20 K and 
another warmer component whose temperature changes from 
25 to 40 K source to source. The cold temperature component 
corresponds to the well-known cold dust in the GC (Odenwald 
& Fazio 1984, Cox & Laureijs 1989, Pierce-Price et al. 2000). 
On the contrary, the detection of a new warm dust component 
with a temperature up to 40 K is a goal of the ISO/LWS spectra 
(see also Lis & Menten 1998, Rodriguez-Fernandez et al. 2000, 
Lis et al. 2001, Goicoechea et al. 2004). Table |4] also shows 
the far-infrared flux of each dust component. It is noteworthy 
that, contrary to the total opacity, the total dust luminosity is 
dominated by the warm dust contribution (90 %; Table|4j. 
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Fig. 4. Sample of LWS01 spectra of selected lines 



2 + 



Fig. 5. LWS04 spectra. For each source the lower spectrum is 
the [Oi] 145 jum line, that in the middle the [Oi] 63 yumline and 
the upper one the [Cn] 158 /mi line. In addition the upper pan- 
els also show 13 CO spectra (gray lines). All the LWS04 spectra 
are represented in units of 10 17 W cm^juirr 1 . The 13 CO spectra 
are represented in arbitrary units. 



To get further insight on the distribution and properties of 
the warm dust in the GC we have compared our results with 
the E band (18.2-25. 1 yum) image of the GC taken by the MSX 
satellite (see Fig. 1). This image shows extended emission over 
the central 300 pc of the Galaxy in addition to a number of point 
sources. We have averaged the MSX E band intensity at the lo- 
cation of our sources in an area comparable to the LWS beam. 
In Fig. we have plotted the temperature of the warm dust 
component derived from our fits versus the MSX E band inten- 
sity. The plot shows a good correlation between both quanti- 
ties, which imply that the MSX image is showing the distribu- 
tion of dust with a temperature higher than 30 K. Apart from 
the well-known radio continuum sources (Sgr A, B, C...), the 
warmest dust is found in the Thermal Filaments and the Sickle 
nebula (both in the Radio Arc region), in a source located at 
(/, b) ~ (0.3, -0.01) and towards (/, b) ~ (-0.35, -0.03). 



4. Fine structure lines 

4.1. Comparison with other observations 

The high spectral resolution of the Fabry-Perot observations 
allows us to study the line profiles and peak velocities of the 
fine-structure lines, especially those of the [Cn] lines, which 
have a high signal-to-noise ratio. To compare the [Cn] data with 
molecular lines observations, in Fig. [5] we overplot the [Cn] 
lines on the I3 CO(1-0) spectra of Rodriguez-Fernandez et al. 
(2001a). In addition, in Fig. [8] we show the central velocities 
of the [Cn] lines (solid triangles) and the linewidths (errorbars) 
derived from the Gaussian fits (Tabled) overlaid in the 12 CO(l- 
0) map of Bally et al. (1987). 

First, it is noteworthy that the [Cn] emission in M- 
0.96+0.13 arises in two components with velocities of and 
-100 km s _1 but no emission is detected from the prominent 
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Table 3. Detailed results of the fits to the dust emission in two 
sources with two independent gray bodies (see text) 



M + 0.21-0.12 



T 


T30tv 




T c 


TiOc 




? 

X~ 


M+0.76-0.05 


25 


0.27 


1 


13.9 


101 


3 


162 


25 


0.27 


1 


14 


16.8 


2 


158 


25 


0.27 


1 


15 


1.83 


1 


157 


24 


0.60 


1.5 


15 


6.3 


1.5 


166 


24 


0.63 
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Table 4. Total flux (FIR) and temperature (T) of the warm (w) 
and cold (c) components needed to fit the spectra with two in- 
dependent gray bodies. 
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1.1 
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33-39 
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8.0 


1.8 


14-15 
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6.6 


0.6 
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26-28 


M+0.76-0.05 
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1.1 


14-16 


23-25 


M+0.83-0.10 


5.5 
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15-16 


25-28 
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0.7 
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0.7 


1.0 


16-19 


27-31 


M+3.06+0.34 


0.3 


0.8 
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31-36 



CO component at 150 km s _1 . In the low spectral resolution H2 
spectra of Rodriguez-Fernandez et al. (2001a) we also found 
evidences that the warm H2 emission arises mainly from the 
negative velocity gas. Thus, the warm H2 and the [Cn] emis- 
sion are probably associated. The sources in the Sgr C area (M- 
0.42+0.01, M-0.50-0.03 and M-0.55-0.05) show clearly two 
[Cn] velocity components. This is in agreement with the hy- 
drogen recombination lines observations of Pauls & Mezger 
(1975). The [Nil] and [Om] data of Rodriguez-Fernandez and 
Martin-Pintado (2004) also show both velocity components. 
However, in contrast to the [Cn], where both components have 

















.X-.> x 



M + 0.76-0.05 




Fig. 6. Dust continuum emission toward M+0.21-0.12 and 
M+0.76-0.05. The black solid points represent the observed 
spectra. For each source we show two different fits with two 
independent gray bodies. For each fit we show the total emis- 
sion in black lines and the emission due to the cold and warm 
components in gray lines. For M+0.21-0.12 the dotted lines 
represent the fit obtained with = 1 and fi c = 2 and the 
dot-dashed lines that obtained with /3 W = 2 and /3 C — 3. For 
M+0.76-0.05 the dotted lines represent the fit obtained with 
p w = 1 and f3 c = 1 and the dashed lines that obtained with 
p w = 1.5 and f3 c = 2.5. The other parameters of the fits can be 
found in Table|5] 



a similar intensity, the [Nn] and [Om] lines are much more 
prominent in the v=-100 km s _1 component, which is gas 
associated to the Sgr C Hn regions. 

In the Radio Arc region (M+0.21-0.12 and M+0.16-0.10) 
the [Cn] lines are centered at ~ 10 km s _1 . Thus, these 
seem to arise in the same gas that was detected by Pauls & 
Mezger (1980) in hydrogen recombination lines. That velocity 
is also in agreement with the [Om] and [Nn] data of Rodriguez- 
Fernandez and Martin-Pintado (2004). On the contrary, the CO 
emission arises mainly at 60 km s _I (Figs.|5]and[S}. Therefore, 
the bulk of the [Cn] emission in the Radio Arc is not associated 
to the dense molecular clouds. Unfortunately, the low spec- 
tral resolution H2 data of Rodriguez-Fernandez et al. (2001a) 
do not allow us to discriminate the origin of the warm H2 
emission. Nevertheless, the warm NH3 emission observed by 
Hiittemeister et al. (1993) in the only source of their sample lo- 
cated in the Radio Arc area (M+0.10-0.12) arises in gas with a 
velocity of 26 km s , between that of the ionized gas and the 
60 kms 1 CO cloud. 
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MSX band E intensity 

Fig. 7. Temperature (K) of the warm grey body needed to fit the 
dust spectra versus the MSX band E (18.2-25.1 fim) intensity 
toward the sources presented in this paper 



For the rest of the sources the agreement between peak 
velocities and line profiles of the [Cn] and CO lines is 
quite good (see for instance M+0.83-0.10 or M+2. 99-0.06). 
However, for some sources the fine structure lines peak veloci- 
ties show a small shift respect to the CO lines (see for instance 
M+0.48+0.03 and M+0.76-0.05 in Fig.0. The magnitude of 
this shift is similar to the Fabry-Perot spectral resolution (30 
km s -1 ) but larger than the wavelength calibration uncertain- 
ties (7 km s _1 ), therefore we believe it is real. The velocity 
shift is also obvious in Fig. [8] where the [Cn] peak velocities 
are always in the envelopes of the dense molecular clouds. This 
finding suggests that the clouds are illuminated from the out- 
side and that the [Cn] emission arises in PDRs in the external 
layers of the clouds. This is in agreement with the recent re- 
sults from the AST/RO survey by Martin et al. (2004). Their 
longitude-velocity diagrams show that the kinetic temperatures 
of the cloud envelopes are higher than those of the central parts 
of the clouds. 

4.2. Extinction 

The fine structure lines presented in this paper have been cor- 
rected for extinction following the method discussed in detail 
in Rodriguez-Fernandez and Martin-Pintado (2004). Basically, 
a lower and an upper limit to the extinction have been derived 
for every source. The lower limits were obtained from the [S in] 
18/33 fim lines ratios and are similar to the foreground extinc- 
tion to the Galactic center as derived from star counts (Ay ~ 25 
mag). On the other hand, conservative upper limits to the ex- 
tinction have been obtained by comparing radio and infrared 
hydrogen recombination lines as well as from the total column 
density of molecular gas measured toward each source. These 
limits range from ~ 45 to ~ 90 mag (Ay) from source to source 
(see Rodriguez-Fernandez and Martin-Pintado 2004). All the 
following analysis has been done taking this range of possible 
extinctions into account. The extinction correction uncertain- 
ties are the main source of errors in the results. To obtain the 
actual values of the extinction at the wavelength of a given line, 



we have used the Draine (1989) law in the mid infrared range. 
For the lines in the far-infrared range, we used the 30 yum ex- 
tinction from Draine (1989). We have extrapolated to longer 
wavelengths using a power law with exponent 1.5. 

4.3. Excitation: shocks versus radiation 

The lines presented in this paper are among the main coolants 
of the interstellar gas with a temperature of a few hundred K. 
Thus, the study of this line emission (in addition to the dust 
continuum of the previous sections) can give us important esti- 
mates of the thermal balance of the GC clouds and the heating 
mechanisms of the warm gas. In the following, we perform a 
similar analysis to that done by Van den Ancker (1999) to study 
the ISO spectroscopy of Young Stellar Objects (Fig[9]has been 
made with data kindly provided by Mario Van den Ancker). 

4.3.1. Shocks 

The left panels of Fig. [9] shows the J-shock model predictions 
(Hollenbach & McKee 1989) for the intensity of several lines 
as a function of the total intensity emitted by all lines for dif- 
ferent shock velocities (30-150 km s _1 ) and preshock gas den- 
sities (10 3 -10 6 cm -3 ). J-shocks models fail to predict the ob- 
served line intensities: first, the observed [Oi] 63 fim intensity 
is considerably lower than that predicted by the models, and 
second, J-shocks models do not predict enough emission from 
[Cn] 158 fim and [Si n] 35 fim even for shock velocities as 
high as 150 km s _1 . 

The center panels of Fig.[9]present C-shock model predic- 
tions (Draine et al. 1983) for the intensity of several lines as 
a function of the total intensity emitted by all lines for shocks 
with velocities of 5-40 km s _1 and preshock gas densities of 
10 4 -10 6 cm 4 . The [Oi] 63 fim intensities can be explained 
by a C-shock with a preshock density in the range of 10 3 -10 5 
cm -3 and shock velocities of 10-20 km s _1 . The detected [Oi] 
146 fim intensities as well as the upper limits suggest densi- 
ties somewhat higher (10 4 - 10 5 cirT 3 ) than those estimated 
with the [Oi] 63 fim line and slightly lower velocities (5-10 
km s _1 ). The [S i] lines seem to rule out C-shocks in gas with 
a pre-shock densitiy higher than 10 5 cm" 3 . In summary, fine 
structure lines from neutral species can arise in C-shocks with 
10 kms 4 and 10 4 cm 4 . 



4.3.2. Diffuse ionized gas and PDRs 

The right panels of Fig.|5]show the intensities of several lines 
versus the incident far-ultraviolet flux as predicted by PDR 
models (Tielens & Hollenbach 1985) with gas densities in the 
range of 10 2 -10 6 cm -3 . As discussed by Van den Ancker 
(1999), in PDRs the far-infrared intensity (Ifir) is correlated 
to the far-ultraviolet incident intensity (Ifuv) by the relation 



'FUV 



0.5 



Q 2 



(3) 



where Q is the spatial extent of the PDR. Here we have as- 
sumed that the PDRs fill the LWS beam. To obtain the intensity 



8 



N. J. Rodriguez-Fernandez et al.: ISO observations of the Galactic center Interstellar Medium: 



in terms of Go, we have divided Ifuv by the average local in- 
terstellar FUV field (1.2 10 4 ergs 1 ctrT 2 sr 1 , Habing 1968). 

The [Fe n] 26 fim and [S i] 25 //m intensities and upper 
limits are compatible with the predictions of the PDR models. 
The intensity of the [Oi] lines can be explained with a density 
of 10 3 -10 4 cm 3 and a FUV incident field G of ~ 10 3 3 
for all the sources except those located in the extremes of 
the Central Molecular Zone (CMZ, Morris & Serabyn 1996) 
and those located in the Clump 2 (M-0.96+0.13, M+ 1.56-0.30, 
M+3. 06+0.34 and M+2.99-0.06). The line emission in these 
sources would arise in PDRs of a similar density (10 3 cm -3 ) 
but with a lower Go (~ 10 2 5 ). That is, the interstellar radiation 
field would be around 10 times lower than in the internal re- 
gions of the CMZ (from Sgr C to Sgr B2). The intensity of the 
[Cn] 158 fim line can also be naturally explained in the context 
of a 10 2 -10 5 cm PDR. However, the PDR model prediction 
for the [Si n] 35 //m line intensity is ~ 10 times lower than 
observed. 

The ionization potential of Ci and Sii is lower than 13.6 eV. 
Therefore, [Cn] and [Si n] lines can arise from mainly neutral 
as well as from ionized regions. The photo-ionized gas com- 
ponent presented in Rodriguez -Fernandez and Martin-Pintado 
(2004) can also contribute to the flux of these lines. To inves- 
tigate this possibility in Fig. ^3 we have plotted the [Cn] line 
flux versus the [Nn] flux (both lines have been observed with 
the LWS beam) and the [Si n] line versus the [S m] 33 yum line 
(both observed with the same SWS beam). The two plots show 
remarkable correlations. Contrary to the [Si n] and [Cn] lines, 
the [Nil] and [S m] can only arise in the ionized gas (the ioniza- 
tion potential of N and S + is higher than 13.6 eV). Therefore, 
the observed correlations imply that a fraction of the [Cn] and 
[Si n] fluxes should indeed come from the ionized gas. 

A least squares fit to the [Cn] and [Nn] data gives the fol- 
lowing result: 

F(CII) = (3.8 + 0.6)F(NII) + (4 ± 2)10~ 18 (4) 

where F(CII) and F(NII) are the [Cn] 158 and [Nn] 122 
yum fluxes in W ctrT 2 . For comparison, Goicoecha et al. (2004) 
found an slope of 5.2 in the Sgr B2 envelope and Malhotra 
et al. (2001) found an slope of 4.3 in their sample of normal 
galaxies. Since the bulk of the [Nn] emission arises in the ex- 
tended low density ionized gas, we conclude that a fraction of 
the [Cn] emission also arises in this gas component. One can 
approximate the [Cn] flux arising from the PDR, F(CIIpdr), 
as F{CU PDR ) = F(CII) - 3.8 F(NII) which is on average of 
4 10 W ctrT 2 . Therefore, for the Clump 2 and the outer 
CMZ sources the [Cn] flux is likely to be mostly of PDR origin 
while for the inner CMZ sources up to ~ 70 % of the flux could 
arise in extended low density ionized gas. 

On the other hand, a least squares fit to the [Si n] and [S m] 
fluxes gives 

F(SiII) = (1.11 ±0.0.09)F(5///) + (4 + 3)10" 19 (5) 

where F{SiII) and F{SIII) are the [Sin] 35 jum and 
[S m] 33 jum fluxes in W ctrT 2 . To understand this rela- 
tion we have computed some CLOUDY (Ferland 1996) mod- 
els. We have studied the [Si n]/[S in] ratio in a range of elec- 
tron densities (0.1-10 4 ctrT 3 ), effective temperatures of the 



ionizing radiation (25000-40000 K) and ionization parameters 
{logU = -1,-4). In that range of parameters, the [Si n]/[S in] 
ratio is almost independent of the electron density and only 
weakly dependent on the radiation temperature. The predicted 
[Si n]/[S m] ratio in the ionized gas is ~ 0.2 - 1 for ionization 
parameters between logU = -2 and -3, respectively, which is 
the range derived from the data of Rodriguez-Fernandez and 
Martin-Pintado (2004). Since the observed correlation gives a 
[Si n]/[S m] ratio in the range 1.0-1.2 one concludes that most 
of the [Si n] emission could indeed arise in the ionized gas. This 
result explains naturally the observed [Si n] excess when com- 
pared to the PDR models. Alternatively, the silicon gas phase 
abundance in the GC could be much higher than that consid- 
ered by the PDR models (~ 3.6 10 ). However, to explain the 
observed flux excess all the silicon should be in gas phase. We 
conclude that a fraction of the [Si n] flux comes from the ion- 
ized gas. For NGC 7023, this result is in agreement with the 
fact that the [Si n] line is more intense toward the ionized gas 
than toward the PDR itself (Fuente et al. 2000). 




[N II] 122 (W/cm2) [S III] 33 (W/cm2) 



Fig. 10. Left: [Cn] 158 //m versus [Nn] 122 jjm line fluxes. 
Right: [Si n] 34 /jm versus [S m] 33 yum line fluxes 

In summary, the fine structure lines of neutral species are 
consistent with a C-shock or a PDR. However, the lines from 
ionic species imply a PDR/ionized gas scenario. 

5. H 2 pure rotational lines 

The H2 pure rotational lines toward the GC clouds have been 
discussed by Rodriguez-Fernandez et al. (2000, 2001a). The 
S(0) and S(l) lines trace a gas temperature of ~ 150 - 200 
K while the lines up to the S(5) probe gas temperatures up 
to ~ 500 K. The column density of H2 at ~ 150 K is a few 
10 22 cirT 2 while the column density of gas at 500 K is only 
~ 1% of that at 150 K. On average, the column density of gas 
with a temperature of 150 K is ~ 30% of the total gas column 
density estimated from CO but there are significant differences 
from source to source (the warm gas fraction in M-0.96+0.13 is 
close to 1). The H2 excitation was also discussed by Rodriguez- 
Fernandez et al. (2001a) by comparing the population diagrams 
derived from the observations with the same type of diagrams 
constructed from J- and C-shock and PDR models. The results 
were that the emission from lines arising from more excited 
levels than those connected to the S(l) line are consistent with 
a C-shock of 12 km s 1 and 10 6 cm 3 , a PDR of G = 10 4 
and 10 6 cirT 3 (both comparing with models and observational 
templates as NGC7023) or a J-shock of 50 km s 1 and 10 6 
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cm -3 . In contrast, to explain the S(0) and S(l) intensities sev- 
eral PDRs with lower density and incident field (Go = 10 3 , 
n — 10 3 cn"T 3 ) or several C-shocks with velocities as low as 7 
km s _1 are required. Of course, a combination of these PDRs 
and low velocity shocks is also possible. 

In this section we revisit the H2 excitation in the same con- 
text in which we have studied the fine structure line emission 
in the previous section. Figure ITT1 shows the intensity of the 
S(0) and S(l) lines as well as the excitation temperatures de- 
rived from the J=3 and J=2 levels (7^2) and the J=3 and J=7 
levels (r 7 3). These quantities are represented versus the total 
intensity radiated by lines and compared to different J-shocks 
(left panel) and C-shocks (middle panel) model predictions. 
The right panel of Fig.[^presents our S(0), S(l), 7/32 and 7/73 
values versus Go (derived from the far-infrared flux). We also 
show a comparison with PDR model predictions. 

The H2 emission is difficult to study in J-shocks due 
to the dissociation and subsequent molecule reformation. 
Furthermore there are no model predictions for the low-J H2 
pure rotational lines in shocks with preshock density lower than 
10 4 cm 3 . The observed 773 and 7^2 are in the range of pos- 
sible values expected by J-shock excitation. However, the total 
intensity radiated by lines in a J-shock with preshock density 
higher than 10 4 ctrT 3 is much higher than that observed from 
the GC sources. In addition, the measured line intensities of 
the S(0) and S(l) lines would be at least 2 orders of magni- 
tude higher than those predicted by these models. In summary, 
J-shocks models cannot explain the observed H2 emission and 
the total energy radiated by lines. 

The central panels of Fig.[^show that both 7/73 and T12 can 
be explained by C-shocks models with densities of 10 4 - 10 6 
cm -3 and velocities of 5-15 km s _1 . Moreover, the total inten- 
sity radiated by lines in such shocks is in agreement with the 
observations. Figure [^ a l so shows that the observed S(l) in- 
tensities are around 3 times higher than those predicted by a 10 4 
cm" 3 and v = 10 km s _1 shock. This effect is much stronger 
for the S(0) line since the observed intensities are around a fac- 
tor of 10 higher than those predicted by the models. 

The right panels of Fig. ITT1 show that both 773 and 732 can 
be explained in PDRs of n = 10 3 — 10 s cm 4 . Moreover the 
observed far-infrared flux due to the warm dust (and Go) is also 
consistent with this picture. However, the observed S(l) line 
intensity is a factor of ~ 10 higher than predicted by a PDR 
model with Go = 10 3 and n — 10 6 cm 3 . The discrepancy 
between observations and models is even higher for the S(0) 
line. 

In summary, the H2 emission from excited levels can be 
explained either by a C-shock or a PDR. However, neither C- 
shocks nor PDR models can account for the large observed in- 
tensities of the S(0) and S(l) lines. 

6. Discussion: the origin of the heating and the 
chemistry of the GC clouds 

6.1. Photo-dissociation regions 

In the previous sections we have shown that the GC clouds ex- 
hibit a warm dust component with a temperature up to 45 K 



that sould be heated by UV radiation. Regarding the gas, on the 
one hand, the fine structure lines of neutral species and the H2 
emission from excited levels could be heated by a C-shock or 
a PDR. On the other hand, the fine structure lines from ionized 
species are best explained in a PDR/DIG (diffuse ionized gas) 
scenario. The large scale observations of the H2 v = 1 - S(l) 
ro-vibrational line at 2.12 //m by Paket al. (1996) has also been 
interpreted as arising in a Go = 10 3 and n = 10 3 cm -3 PDR. 
These parameters are in excellent agreement with those derived 
in the previous sections. Clearly, the extended PDR/DIG sce- 
nario is the easiest way to explain all the observational data. 
Namely, the warm dust, the fine structure lines from neutral 
and ionized gas, and the H2 emission from excited levels. 

As discussed in Rodriguez-Fernandez and Martin-Pintado 
(2004), the interstellar medium in the 300 central pc of the 
Galaxy is quite non-homogeneous and it is permeated by rel- 
atively hot stellar radiation (~ 35000 K). The DIG emission 
seems to arise from the external layers of the dense molecu- 
lar clouds. Indeed, molecular, neutral/atomic and ionized gas 
are associated (Sect. 4.1) and therefore a PDR is expected in 
the interface between the dense molecular gas and the ionized 
gas. In the external layers of such a PDR, the gas is heated 
by photoelectric effect on the dust grains to temperatures well 
upon the dust temperature. For instance the standard model of 
a Go = 10 3 and n = 10 3 cm -3 PDR by Hollenbach et al. 
(1991) shows that in the first two or three magnitudes of visual 
extinction inside the PDR, the gas is heated to a temperature 
of 100-200 K while the dust temperature is 30-50 K. Thus, the 
discrepancy of gas and dust temperatures is normal in PDRs. 
Furthermore, the measured gas and dust temperatures in the 
GC clouds are in excellent agreement with the Hollenbach et 
al. (1991) models. However, the column density of molecular 
gas with a temperature of ~ 150 K predicted by these mod- 
els is 2-3 10 21 cirT 2 , i.e., 10-30 % of the measured warm gas 
column density in the GC clouds (Rodriguez-Fernandez et al. 
2001a). For some lines of sight the discrepancy could be lower 
if one takes into account that several clouds are observed (dif- 
ferent velocity components) and each of them can present a 
PDR. Nevertheless, even with 2 or 3 standard PDRs in the line 
of sight, it will be difficult to explain the large column densities 
of warm gas measured in the GC clouds. 

In addition, Rodriguez-Fernandez et al. (2001a) have mea- 
sured a warm NH3 abundance in the GC clouds of ~ 10~ 7 , 
which is only 4 times lower than the typical C lfi O abundance 
in the GC (Wilson & Matteucci 1992). This large NH 3 abun- 
dance is difficult to explain in the context of a PDR. For in- 
stance, the NH3 to C 18 abundance ratio measured in the Orion 
bar PDR is ~ 10~ 2 (Batrla & Wilson 2003), well below the 
value measured in the GC. In a PDR the warm NH3 emission 
arises from the external layers, where the gas is more exposed 
to the radiation. The problem is that the NH3 molecule is easily 
photo-dissociated by UV radiation. The NH3 dissociation rate 
can be expressed as ~ 3 10 _10 G exp[-2.16 A v + 1.71 10~ 2 A 2 ] 
s _1 (Stenberg & Dalgarno 1995). Therefore the dissociation 
rate in the warm region (Ay < 2 mag.) of a Go = 10 3 PDR is 
~ 4.3 10~ 9 s . Thus, there would be no warm NH3 after ~ 7 
yr. One can estimate whether the NH3 abundances could be ex- 
plained in a dynamical PDR in which NH3 rich gas is continu- 
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ously been injected in the photo-dissociation front. To maintain 
a constant NH3 abundance of 10~ 7 , a NH3 column density of 
2 10 14 cirT 2 (H2 column density of 2 10 21 ctrT 2 ) should be 
injected to the PDR every 7 yr. Taking into account a H 2 den- 
sity of 10 3 ctrT 3 one derives that the NH 3 rich gas should be 
entering the photo-dissociation front with a velocity of several 
thousand km s _1 , which is not observed. 

The large warm H2 column density and the NH3 abundance 
measured in GC clouds is difficult to explain in a PDR scenario, 
even if the photo-dissociation front is not stationary. In the fol- 
lowing we review other alternatives. 

6.2. Cosmic rays 

A cosmic rays ionization rate of ~ 10~ 17 s _1 maintains the 
gas temperature at ~ 10 K in the local quiescent molecular 
clouds. An ionization rate 100 times higher than the local value 
would balance the gas cooling rate at 70 K. Indeed, an en- 
hanced cosmic rays ionization rate was one of the first mech- 
anisms proposed to explain the warm temperatures of the GC 
clouds (Giisten et al. 1981, Morris et al. 1983). However, it is 
now clear that the GC clouds present a wide rage of tempera- 
tures from 20-30 K to 150 and 500 K (Hiittemeister et al. 1993, 
Rodriguez-Fernandez et al. 2001a). Temperatures of ~ 70 K, 
which are sometimes cited as typical of the GC clouds, are just 
averaged values without much physical meaning. Heating by 
cosmic rays is a very pervasive mechanism since cosmic rays 
penetrate and traverse the clouds. If the cosmic rays ioniza- 
tion rate is high enough to heat the gas to 150 K most of the 
gas in the GC clouds should present this temperature. On the 
contrary, Rodriguez-Fernandez et al. (2001a) showed that only 
30% of the gas is warm. Furthermore, EGRET observations do 
not show any excess of y-rays associated to the GC molecular 
complexes that could be due to a cosmic rate ionization rate 
higher than the local value (Mayer-Hasselwander et al. 1998). 
We conclude that cosmic rays cannot be the heating mechanism 
of the bulk of the 150 K gas. 

6.3. X-rays dominated regions (XDRs) 

The GC is an extended source of Fe 6.4 keV emission (Koyama 
et al. 1996). This line arises from neutral or weakly ionized gas 
illuminated by a hard X-ray continuum. Martin-Pintado et al. 
(2000) have found a spatial correlation of the SiO(l-O) and the 
Fe 6.4 keV emission that suggest that both lines have a related 
origin and therefore, it extends the possible influence of the 
X-rays to the molecular gas. The GC clouds could be X-ray 
Dominated Regions (XDRs) as those modeled by Maloney et 
al. (1996). X-rays act directly by photo-electric effect on the 
gas and their attenuation is directly proportional to the gas col- 
umn. As a result, the column density of warm gas in XDRs 
can be 10 times higher than that of a PDR (in PDRs the domi- 
nant gas heating mechanism is the photo-electric effect on the 
dust grains and the UV absorption by the dust is exponential). 
Therefore, XDRs could heat large quantities of gas, such as 
those measured in the GC. 



XDRs can be characterized by the effective ionization pa- 
rameter (e e ff) and the gas density (Maloney et al. 1996). In 
clouds with a density of 10 3 ctrr 3 and a column density of 
10 22 ctrT 2 , e e ff can be expressed as: e e fj = 1.1 \0 6 L^r^ 2 , 
where L 34 is the X-rays luminosity in units of 10 34 erg s _1 and 
r is the distance from the cloud to the continuum source in pc. 
The X-rays luminosity in scales of 8' (~ 20 pc) around SgrB2, 
C or D is ~ 0.3 10 35 erg s _1 while around Sgr A it is approx- 
imately 3 times higher (Sidoli et al. 1999). Taking L34 = 100 
and a spatial scale of 1 pc one gets e e fj = 1.1 1 0~ 6 . In contrast, 
e e //;i0~ 3 is required to heat the gas to temperatures higher than 
100 K (Maloney et al. 1996). Therefore, XDRs do not seem to 
be a possible heating mechanisms in large scales for the bulk 
of the gas. 

6.4. Turbulence or low velocity shocks 

The sound velocity in a gas with a temperature of ~ 100 K is 
~ 1 km s _1 . Hence, the large widths of the molecular lines ob- 
served towards the GC (10-15 km s~') imply the presence of 
a high degree of supersonic turbulence. In addition to the large 
line-widths and the widespread high gas temperatures, the GC 
clouds exhibit a particular hot-core-like chemistry. High abun- 
dances and extended emission of molecules like C2H5OH, SiO 
(Martm-Pintado et al. 1997, 2000, 2001), NH 3 (Hiittemeister et 
al. 1993, Rodriguez-Fernandez et al. 2001a), SO, CH 3 OH (Lis 
et al. 2001) are measured towards large regions (> 10 pc) in 
the whole GC. These molecules are usually in the grain man- 
tles. The relatively low dust temperatures in the GC and the 
extended emission of those molecules imply that the processes 
giving rise to the measured high gas-phase abundances are dif- 
ferent than those taking place in a hot-core, where they are ex- 
plained by thermal evaporation from dust grains with a temper- 
ature higher than 70 K. Furthermore, as discussed above, these 
molecules are easily dissociated in the presence of UV radia- 
tion. Therefore, their high gas phase abundances and their ex- 
tended distribution suggest that they have been ejected from the 
grains by mechanical processes, such as low velocity shocks. 
These have also been invoked to explain the transitory gas heat- 
ing and the non-equilibrium H2 ortho-to-para ratio measured in 
two GC clouds (Rodriguez-Fernandez et al. 2000). 

The mechanical processes could be turbulence or shocks. 
Rodriguez-Fernandez et al. (2001a) showed that the energy in 
turbulent motions is enough to explain the cooling of the 150 
K gas by CO and H2 rotational lines. It is quite difficult to deal 
with the turbulence. Current models do not make direct predic- 
tions of lines intensities but explain the processes taking place 
in the turbulent interstellar medium. For instance, Smith et al. 
(2000) have shown how the turbulence creates a full spectrum 
of shocks with different velocities. They found that the number 
of shocks is inversely proportional to the square of the shock 
velocity. Hence, most of the gas mass is suffering slightly su- 
personic shocks with low Mach numbers of 1-5, which for the 
GC warm gas implies very low shock velocities of 1-5 km s _1 . 
Indeed, the velocity of the shocks taking place in the GC must 
be low enough to heat the H2, but without giving rise to much 
fine-structure line emission since, as we have already shown, 
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the observed fine-structure emission can be accounted for by 
the PDR component alone. The key in the GC could be the 
magnetic field strength, which can be as high as 1 mG. In the 
presence of a strong field C-shocks have an important magnetic 
precursor that heats large quantities of gas to moderate temper- 
atures (~ 150 K). For instance, a C-shock with velocity of 25 
km s _1 in a preshock gas density of 10 4 cirT 3 can reproduce 
the S(0) and S(l) H 2 intensities observed in the GC without giv- 
ing much emission in fine structure lines ([Oi] and [Cn] lines 
intensities lower than 10~ 6 ergs -1 cm _2 sr _1 ; S. Cabrit priv. 
communication). 

The origin of the turbulence and shocks in the GC can be 
local or related to the large-scale properties of the Galaxy. On 
the one hand, massive stars are an important source of shocks 
and turbulence in scales of a few pc by means of supernova 
blast waves, Hn regions and strong winds. For instance, hot 
(60-130 K) NH 3 shells with sizes of 1-2.5 pc have been ob- 
served in the envelope of the Sgr B2 molecular cloud (Martin- 
Pintado et al. 1999). These shells expand with a low veloc- 
ity of ~ 5 km s _1 and are probably created by the interac- 
tion of Wolf-Rayet stars winds with the surrounding ISM. The 
scenario found in the Sgr B2 envelope would explain the non- 
homogeneity of the GC ISM required to explain the gas ioniza- 
tion (Rodriguez-Fernandez et al. 2001b, Rodriguez-Fernandez 
and Martin-Pintado 2004). In addition these findings have im- 
portant implications in the heating of the GC molecular gas. 
The UV radiation from the stars responsible for the ionization 
will also create PDRs in the edges of the cavities. In addition, 
C-shocks with a velocity similar to the shell expansion are also 
expected. 

On the other hand, the origin of the turbulence can be tidal 
forces induced by the rotation of the Galaxy, as proposed by 
Wilson et al. (1982) to account for the high temperatures of 
the molecular gas in the Sgr B2 envelope. Indeed, the models 
by Das & Jog (1995) have shown that tidal forces can heat the 
molecular clouds moving along elongated orbits in a bar po- 
tential, which seems to be the case in the GC (Binney et al. 
1991). 

The GC is a complex region, where many processes can 
play a role in the excitation of the ISM. We conclude that, ex- 
tended PDRs can explain the fine structure line emission and 
the H2 emission from excited levels. On the other hand, the 
most likely heating source of the bulk of 150 K molecular gas, 
traced by the H2 emission from low levels (and the GC chem- 
istry) are moderate velocity shocks most likely induced by the 
turbulent nature of the GC ISM. 



7. Conclusions 

We have presented fine-structure line data emitted by the main 
coolants of the neutral gas and the dust continuum spectra 
around the emission peak. The spectra were taken towards a 
sample of molecular clouds distributed along the 500 central 
pc of the Galaxy and located far from thermal radio continuum 
or FIR sources. We also compared these results with the H2 
pure-rotational lines presented by Rodriguez-Fernandez et al. 
(2000,2001a). 



The dust continuum emission can be modeled with two 
temperatures: a cold component with a temperature of ~ 15 
K and a warmer component whose temperature varies between 
25 to 40 K from source to source. 

We have compared the fine-structure lines and the dust con- 
tinuum emission with the predictions from shocks and PDR 
models. We have found that the observations are best explained 
by a PDR with a density of 10 3 cm 3 and a far-ultraviolet in- 
cident field 10 3 times higher than in the local ISM. The PDRs 
should be located in the interfaces between the dense molecular 
matter and the ionized gas component discussed by Rodriguez- 
Fernandez and Martin-Pintado (2004). Indeed, the [Cn] 158 
yum and the [Si 11] 34 /im lines have a significant contribution 
from the ionized gas component. 

The H2 pure-rotational emission from excited levels 
(Rodriguez-Fernandez et al. 2000, 2001a) and the ro- 
vibrational near-infrared emission (Pak et al. 1996) also arise in 
the PDR component. However, only 10-30 % of the gas traced 
by the lowest pure-rotational levels (S(0) and S(l)), with a tem- 
perature of ~ 150 K, can be accounted for in this context. We 
suggest that the most promising heating mechanism of this gas 
(and the origin of the chemistry characteristic of the GC) are 
low velocity C-shocks created by turbulent motions. 

Higher spatial and spectral resolution as well as mapping 
observations of the warm H2 and key fine-structure lines from 
neutral and ionized atoms are needed to get further insight in 
the interplay of the different ISM phases and the heating of 
the molecular gas in the GC. The new generation of medium 
and far-infrared telescopes (VLT, SOFIA, Herschel, ALMA) is 
well suited to attack these problems in the GC and other nearby 
galaxies. 
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Fig. 8. Peak velocities (triangles) and linewidths (errorbars) of the [Cn] 158 fim line overplotted on the CO(l-O) longitude- 
velocity diagram by Bally et al. (1987) 
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Fig. 9. Left panels: Fine-structure lines intensities versus the total intensity radiated by lines. Solid and empty squares represent 
the observed results for the low and high extinction corrections, respectively. For comparison, the lines intensities versus the 
total intensity radiated by lines in J-shocks of different velocities (dashed lines) and densities (dotted lines) are also shown. 
Central panels: Same for C-shocks of different velocities (dashed lines) and densities (dotted lines). Right panels: Intensities 
of different lines versus the far-ultraviolet incident flux. Solid and empty squares represent the results for the low and high 
extinction corrections, respectively. For comparison, the PDR model predictions for different densities are also shown 
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Fig. 11. Left panel: H2S(0) and HaSfl) intensities and rotational temperatures derived from the J=3 and J=2 levels (T32) and from 
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